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Abstract. The time evolution behaviour of the superhumps of the dwarf nova 1RXS J232953. 9+062814 is in- 
vestigated with the wavelet analysis method. On the basis of two nights CCD photometry performed during its 
first superoutburst as well as other published brightness data, we reveal the superhump's time-dependence as a 
function of periods and time. Our light curves, which phased in the rapid decay ending portion of the superout- 
burst and in the dawn of a following normal outburst, are important to help trace the superhump evolution for 
the star. Evident amplitude variations of the superhumps, reflecting the fading of outbursts, are detected. The 
general profile of brightness fading over the outbursts roughly followed an exponential decay law or a form of a 
five-order polynomial. Both the superhump period and the orbital period of the binary system are detected in 
the present data. We obtain P sh =0.04575±0.00005 d and P orb =0.04496±0.00005 d. They agree with the existing 
values based on additional data. The two periods exchanged their roles during the superhump evolution. 
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1. Introduction 

Cataclysmic variables (CVs) are close binary systems in 
which a low-mass secondary transfers mass onto a white 



dwarf (Warner 1995), while dwarf novae are a sub-group 
of CVs which experience repetitive outbursts with typi- 



cal amplitudes of 2-5 mag (Warner 1987). The dwarf nova 
outbursts are considered to be a suddenly enhanced re- 
lease of gravitational energy induced by the thermal in- 
stability of an accretion dick (Osaki 1996| ). During out- 
bursts, short-period dwarf nova systems often show pho- 
tometric oscillations (superhumps) at periods of a few per- 
cent longer than their orbital periods (P rb)- The super- 
hump period is thought to be the beat between P or b and a 
tidally driven precession of the eccentric disk. According 
to the standard evolut ion mo del o f compact binary sys- 
tems (Paczynski 1981 ; King 1988 ), orbital angular mo- 
mentum losses maintain the mass-transfer and thus drive 
the CV evolution. 

1RXS J232953.9+062814 (equinox at 2001.11.03.926) 
was i dentified as a cataclysmic variable by Wei et al. 
( 1999 ). It was classified as a dwarf nova by Hu et al. 
( 1998 ). Its outburst was first detected on 2001 November 
3.926 UT by P. Schmeer, a member of the VSNET 
collaboration teamQ This team soon reported super- 
humps with amplitudes of 0.2-0.3 mag and a period of 



0. 046311+0. 000012d, indicating that this object is a dwar f 
nova of type SU Ursae Majoris (Uemura et al. 2001 ). 
I n ter ms of the quiesce nt sp ectra obtained by Hu et al. 
( 1998 ) and Wei et al. ( 2001 ), this object is a hydrogen- 
rich CV with a high-inclination accretion disc. However, 
the observed short superhump period is below the common 
'period minimum' for hydrogen-rich secondaries (~1.3h). 
Further data showed that the star has a short distance 



of 140-350 pc (Skillman et al. |2002j; Wei et al. |200l|) and 
it is much brighter than V485 Centauri, the other known 
object of this class. Therefore, the dwarf nova is a very 
important object for studying the evolutionary scenario 
of CVs and it was immediately monitored by two research 
networks after its detection of outburst with the aim to 
detect the time-evolution of the superhumps and the sys- 
tem's properties (Skillman et al. 2002; Thorstensen et al. 
[2002 ; Uemura et al. 2002). Our observations were carried 
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out at the Xinglong Station of the National Astronomical 
Observatories of China (NAOC). This paper reports the 
results of a study on the superhump evolution based on 
all available brightness data. 

2. Photometry 

One night coordinated observations of the dwarf nova 
1RXS J232953.9+062814 were carried out on two tele- 
scopes at the Xinglong Station of NAOC on 2001 
November 8. The Johnson V photometry was performed 
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Fig. 1. The CCD image of the filed of the dwarf nova 
1RXS J232953. 9+062814, obtained with the 60-cm super- 
nova survey telescope at the Xinglong Station of NAOC. 



with the three-channel high-speed photoelectric photome- 
ter, designed fo r the Whole Earth Teles cope campaign 
(Nather et al. |199C| ; Jiang & Hu |1998| ), mounted on 
the 85-cm Cassegrain telescope. The star GSC 0591- 
1706 (a 2 ooo = 23 h 30 m ll s .42, <S 2000 = 06°26'07".8, V = 
10.52 mag, K2); was used as comparison as suggested by 
VSNET members. The dwarf nova, comparison and sky 
background were simultaneously exposed in continuous 5- 
s intervals. The unfiltered CCD photometry was carried 
out on the 60-cm telescope, which dedicated to a super- 
nova survey. Another night CCD photometry was contin- 
ued on 2001 November 9. The CCD camera was made by 
Princeton Instruments, which is a liquid nitrogen cooled 
CCD and it has a peak quantum efficiency of over 90%. It 
has 1340 x 1300 pixels, each 20/im in size. The field of view 
is 10 ' x 10' at the telescope's Cassegrain focus (Qiu et al. 
2001 ). The integration time for each CCD frame is 10 s. 
In total, we obtained 3277 V measurements together with 
1785 CCD frames of the dwarf nova field. Figure |l| shows 
the CCD image of the field. The standard CCD aperture 
photometry was made using the procedures implemented 
in the packages of IRAF. 

Unfiltered CCD differential magnitudes were estab- 
lished with respect to an assemble comparison star con- 
sisting of several stars (Cl^GSC 0591-1660, 14.1 mag; 
C2=GSC 0591-1692, 14.4 mag; C3=GSC 0591-1665, 14.0; 
C4: a 20 oo = 23 h 29 m 42 s .31, S 200 q = 06°26'10".2, USNO, 
14. 5R) in the observed field (Fig j^) . The typical obser- 
vational accuracy for the CCD data is ^0.01 mag. From 
panel 'a' of Fig. ||, the standard deviations for the three 
sets of differential CCD measurements (C1-C2), (C1-C3) 
and (C1-C4) are 0.010, 0.013 and 0.018 mag, respectively. 
The light curves are presented in Figs. || and ||, where the 
fits (solid and dashed lines) with the superhump period 
and orbital period are also plotted. 



The photoelectric data, unfortunately, had some prob- 
lems. It seemed that the sky and comparison channels 
were abnormal during the observations. We failed to use 
them to establish differential magnitudes for the variable. 
Finally, we used the variable's measurements alone after 
subtracting a second-order polynomial fit. The results are 
only valuable for a rough estimate of the light variations 
in the star. The light curves have been given in Wei et al. 
(2001). For this reason, the data were not considered in 
the following analyses. 



3. Wavelet analysis 

In the analysis of time-series data of variable stars, be- 
yond the detection of periodicities, one needs to investi- 
gate the stability of the periodicities as the oscillations 
are caused by matter transfer from a component to the 
other as in CVs. The detection of a variable period is very 
important as it gives information on the binary nature or 
on the evolutionary status of a star. However, it is diffi- 
cult to determine the period variability. The conventional 
Fourier analysis is unable to give any information about 
period variation due to it just localizes frequency infor- 
mation without time localization. Unlike Fourier analysis, 
in which we analyze signals using sines and cosines, alter- 
natively, wavelet analysis uses wavelet basis to decompose 
and reconstruct a periodic signal or an oscillation. The 
wavelet basis localizes both in time and frequency spaces, 
so the wavelet analysis is a well suited method to accom- 
plish the mission of the detection of time dependence of 
oscillation frequencies. In fact, it has been proving to be a 
versatile tool to disclose periodic and time-dependent be- 
haviour of a signal. Moreover, wavelet transforms do not 
ask a long-time coverage data set as Fourier transforms 
do in resolving periodicities in a signal. 

There are a quite number of examples using wavelet 
analysis method in the study of var ious oscillating be- 
haviour s. Fo r instance, Gou pil et al. ( 1991 ), Szatm ary fc 
Vinko ( |1992D , Sca rgle e t al. fll993j ), Coupinot et al. ( |1992 | ) 
and Starck et al. ( |1997|) . In special, Fritz & Bruch ( 1998 ) 
used wavelet transforms to study the flickering light curves 
in CVs. The applicability of wavelet analysis for studying 
periodicity in the light curves of variables has been fur- 
ther explored by some astronomers such as Szatmary et 
al. ( |1994| , |1996| ) and Foster ( |1996| ). The wavelet transform, 
as these authors have demonstrated, is well-suited to de- 
tect the local behaviour of the light curves of variables. For 
example, the investigation of time-dependent phenomena 
including amplitude and frequency modulation, changes 
of period and phase. In view of the unevenly spaced data 
points, the ability of several improved wavelet transforms 
to detect and quantify periodic and pseudo-periodic sig- 
nals, the weighted wavelet Z-transform (WWZ) developed 
by Foster (1996) was used. 

We set trial detection for frequencies from 8 to 50 cy- 
cle d _1 (corresponding to periods in the range of 0.02 
to 0.125 d). Times varied from to 0.2 d in a step of 
0.002 d. Frequency resolution (or calculation error) was 
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Fig. 2. Light curves on 2001 November 8. (a): Differential magnitudes between the comparison stars CI, C2, C3 and 
C4. Means were subtracted. For display, the means of (C1-C3) and (C1-C4) were shifted to 0.08 and 0.16 mag from 
zero, respectively, (b): Differential magnitudes of 1RXS J232953. 9+062814 relative to CI, (V-Cl), and a linear fit: 
0.924 ± 0.005 + (2.125 ± 0.044) t,t =JD~2452221.9667d. (c): (V-Cl) but the linear trend was removed, along with 
sinusoids of periods P s h (plus 0.5P s h, dashed line) and of both P s h and P or b (solid line). 



1.0 cycle d 1 . This tended to establish rough limits for 
the fluctuation frequency as a function of time. We ob- 
tained frequencies at 23, 36, 37, 38, 39, 40, 41, 22, 44, 43, 
42 and 41 cycle d _1 in order with time increase for the 
2001 Nov. 8 data, at which the WWZ transform reaches 
maxima. As for the 2001 Nov. 9 data, we obtained 41, 
42, 43 and 22 cycle d _1 . Most of the data points at most 
times were of frequency 22 cycle d _1 , a value close to 
what we want to find, the superhump period 0.0463 d 
(21.59 cycle d" 1 ) and the orbital period 0.0446(22.42 cy- 
cle d _1 ), respectively (see the distribution in Fig. ||). At 
the first place, we see the frequency or period was not 
constant during the observing span. In addition, the two 
sets of trial periods appeared in opposite phases. This is 



why we use wavelet analysis rather than Fourier analysis. 
Fourier transform always assumes periods to be constant 
in the investigated duration, so it prevents us from un- 
covering the variability of a period. Then we refined the 
step-length of frequency to be 0.02 cycle d -1 , correspond- 
ing to an error of ^0.00005 d in period for the star, to 
obtain more accurate values of periods. We inferred that 
the orbital period (0.0446 d, frequency 22.42 cycle d _1 ) or 
its daily alias (21.42) and harmonic (44.84) were mainly 
present in the first night data, while the superhump pe- 
riod (0.0463, frequency 21.60 cycle d _1 ) or its daily alias 
(22.60) and harmonic (43.20) were mainly present in the 
second night data. No trace of 22.42 (orbital frequency) 
was found in the second data set. The second night data 
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Fig. 3. Light curves on 2001 November 9. (a): Differential magnitudes between the comparison stars CI, C2, C3 and 
C4. Means were subtracted. For display, the means of (C1-C3) and (C1-C4) were shifted to 0.08 and 0.16 mag from 
zero, respectively, (b): Differential magnitudes of RX 2329+06 relative to CI, (V-Cl), and a linear fit: 1.585 ± 0.006 — 
(0.520 ±0.043) t,t =JD-2452222.9412d. (c): (V-Cl) but the linear trend was removed, along with sinusoids of periods 
P s h (plus 0.5P s h, dashed line) and of both P s h and P or b (solid line). 



mainly reflect the superhump content oscillating at 21.59. 
Finally, according to the distribution of periods and num- 
bers of data points where WWZ values reached maxima 
(Fig. H), we obtained frequencies at 22.24±0.02 cycle d _1 
(P orb =0.04496±0.00005d) representing the orbital motion 
of the binary system from the 2001 Nov. 8 data, and at 
21.86 cycle d" 1 (P sh =0.04575±0.00005d), the photomet- 
ric oscillations of the dwarf nova. The results agree with 
the values 0.0446d (64.177 m) and 0.04637(4 ^ (6 6.672 
m), respectivel y, gi ven by Skillman et al. ( 2002 ) and 
Uemura et al. (|2002|) . 

The results of wavelet analyses including the variation 
of superhumps as a function of both time and frequency, 
the time-evolution of amplitudes and periods of the super- 



humps are presented in Figs. | 
curves (see panels 'b' of Figs. 



to H. We used the raw light 
and [}]) without linear-fits 



removed in the analyses. This ensures all information in 
the brightness was not lost. 
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Fig. 4. The distribution of the numbers of effective data 
points at which WWZ reaches maxima and we regard the 
corresponding periods as the main contents of the light 
variations at a time. 



Fig. 6. The wavelet solution of the superoutburst on 2001 
Nov. 9. It shows the time evolution of the superhumps as 
a function of time and frequency (period). 




Fig. 5. The wavelet solution of the superoutburst on 2001 
Nov. 8. It shows the time evolution of the superhumps as 
a function of time and frequency (period) 
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Fig. 7. Time evolution of period, amplitude as deter- 
mined by wavelet transform for the data on 2001 Nov. 



Fig. 8. Time evolution of period, amplitude as deter- 
mined by wavelet transform for the data on 2001 Nov. 
9. 



4. Results and Discussion 

We analysed the superhump evolution in the interacting 
binary system 1RXS J232953.9+062814 based on the two 
nights data (2001 Nov. 08-09) obtained by speed CCD 
photometry. Wavelet analyses reveal that the periodicities 
in the light curves changed with time (see the top panels 
of Figs. and|||) because the exchange of the roles of the 
superhump and orbital motion. Superhump is prominent 
during outbursts while orbital motion of the binary system 
is conspicuous in quiescence. In our case, superhump pe- 
riod is noticeable on 2001 Nov. 9, when the star went into 
a new normal outburst, while the orbital period is remark- 
able on 2001 Nov. 8, when the superoutburst evolved to its 
ending portion of fading. According to t he Fo urier solution 
of the light curves by Skillman et al. ( 2002 ) and the so- 
lution using Phases Dispersion Minimization by Uemura 
et al. ( 2002 ), the periodicities involved in the superout- 
burst are mainly contributed by the superhumps, while in 
quiescence the light variation mainly reflects the orbital 
motion of the binary system. In other words, superhump 
period dominates the superoutburst while orbital period 
contribution governs the quiescence. Therefore, it should 
be that there is a time containing more information on the 
two kinds of variations (orbital and oscillating). This time 
is now clearly between our observations, i.e. between JD 
2452222.0 and 2452223.0. Our observations fell into the 
phase near the tail part of rapid decline of the superout- 
burst and in sight of the first subsequent normal outburst. 



Figure g shows the present data in dots (pointed by ar- 
row s). Th is figure was reproduced from fig. 1 of Skillman et 
al. ( 2002 ), which can also be referred to fig. 1 of Uemura et 
al. ( 2002 ). Furthermore, from the middle panels of Figs. || 
and |3|, it is obvious that the brightness evolved in oppo- 
site directions over the two nights, i.e. on 2001 November 
8 the star faded gradually but on November 9 it turned 
to be brighter and brighter. In addition to the overall pro- 
file of the mean brightness changes, we have noticed the 
linear decay trends involved in the nightly light curves. 
The linear changing contents for these two sets of data: 

0. 924 ± 0.005 + (2.125 ± 0.044)<,i =JD-2452221.9667 d 
and 1.585 ± 0.006 - (0.520 ± 0.043)i, t =JD-2452222.9412 
d clearly show that our observations were in the transition 
phase of the superoutburst. The change rate for these two 
nights are different: the former is faster than the latter. 
The fact is that the star really run into its second brief 
outbursts (or an 'echo outburst') from JD 2452223.0 d, 

1. e. 2001 November 9, UT 12:00 (see fig. 1 of Skillman et 
al. 2002). This indicates the value of the present data for 
tracing the star's outbursts evolution. It also explains the 
two known periods resolved by the wavelet analyses based 
merely on our data. 

Our high-speed photometry shows the detailed light 
variations of the superoutburst of dwarf nova 1RXS 
J232953. 9+062814, such as the linear decay content of 
brightness, the superhumps and the doubled-humped 
wave of orbital motion. Because both the superhump pe- 
riod and orbital period are involved in the light curves, 
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Fig. 9. Mean light curves of 1RXS J232953.9+062814 
during the 2001 outbursts. This is a reproduction of fig. 1 
of Skillman et al. (2002) together with our data (two dots, 
indicated by arrows) and two decay fits. 



either the oscillation or orbital motion cannot reproduce 
the observed light curves well. In Figs. ||and |], the dashed 
curves are the sinusoids of periods P s h=0. 0463d (super- 
hump) and 0.5P s h, solid curves are of periods P s h, 0.5P s h 
plus P orb =0.0446d (orbital). 

Wavelet analyses also reveal the details of the fading of 
the outbursts. The bottom panels of Figs. f|] and || display 
the changes of the amplitudes of corresponding periods. 

We may get an insight into the brightness changes in 
a general point of view. Therefore, we made a first-order 
exponential decay fit and a five-order polynomial fit to all 
the mean light curves available (fig. 1 of Skillman et al. 
2002| ) with the exception of the two normal outbursts. The 
results indicate that the brightness decayed following an 
exponential law of y(mag.) = 16.326 ± 0.064 - (168.34 ± 
60.45)e _t ^ 4 - 693±0 - 408 ' or in a five-order polynomial form of 
y = -30.570±8.706 + (5.110±1.205)t-(0.219±0.064)f 2 + 
(0.00457 ± 0.0016)* 3 - (5 ± 2) x 10~ 5 i 4 + (1.87 ± 0.91) x 
10~ 7 i 5 . Here t= JD-2452200 (d). They are displayed in 
Fig. ^, where the dashed line for polynomial fit while the 
solid line for exponential fit. 
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